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Abstract 

Within the framework of QCD factorization (QCDF), power corrections due to penguin annihilation can 
account for the observed rates of penguin-dominated two-body decays of B mesons and direct CP asymme- 
tries A C p(K~Ti + ), A C p(K*~K + ), A C p(K~p°) and A C p(n + Tl~). However, the predicted direct CP-violating 
effects in QCDF for B — ► K~7t°,K~ri, n~r] and B° — ► 7r 7r° are wrong in signs when confronted with ex- 
periment. We show that subleading 1/mj power corrections to the color-suppressed tree amplitude due to 
spectator scattering or final-state interactions will yield correct signs for aforementioned CP asymmetries 
and accommodate the observed k°k° and p°7r° rates simultaneously. Implications are discussed. 
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1. In the heavy quark limit, hadronic matrix elements can be expressed in terms of certain nonperturba- 
tive input quantities such as light cone distribution amplitudes and transition form factors. Consequently, 
the decay amplitudes of charmless two-body decays of B mesons can be described in terms of decay con- 
stants and form factors. However, the predicted rates for penguin-dominated B — ► PP, VP, VV decays (P and 
V denoting pseudoscalar and vector mesons, respectively) are systematically below the measurements (see 
the second column of Table H for a review, see 111]). 1 Moreover, the calculated direct CP asymmetries for 
B° — > K~n + ,K*~n + , B — > K~p Q and B° — ► 7T + 7T~ are wrong in signs when confronted with experiment as 
shown in the same Table. This implies the necessity of taking into account 1 /nib power correction effects. 
In the QCD factorization (QCDF) approach power corrections often involve endpoint divergences. For 
example, the hard spectator scattering diagram at twist-3 order is power suppressed and posses soft and 
collinear divergences arising from the soft spectator quark and the I /nib annihilation amplitude has end- 
point divergences even at twist-2 level. Since the treatment of endpoint divergences is model dependent, 
subleading power corrections generally can be studied only in a phenomenological way. While the endpoint 
divergence is regulated in the pQCD approach by introducing the parton's transverse momentum yfl, it is 
parameterized in QCD factorization as 

W^ln^lWH (1) 
Jo y At, 

for penguin annihilation contributions with A/, being a typical scale of order 500 MeV. 

In the so-called "S4" scenario of QCDF [4] with some appropriate choice of the parameters p A and <j>A, 
the above-mentioned discrepancies are resolved in the presence of power corrections due to the penguin 
annihilation topology. However, a scrutiny of the QCDF predictions reveals more puzzles in the regard 
of direct CP violation. When power corrections due to penguin annihilation are turned on, the signs of 
Acp in B — ► K~7i°, K r], 71 T] and B° — ► 7i°7i° will also get flipped in such a way that they disagree 
with experiment (see the third column of Table HI). The so-called Kn CP-puzzle is related to the difference 
of CP asymmetries of B — ► K~n° and B° — ► K~n + . This can be illustrated by considering the decay 
amplitudes of B — > Kn in terms of topological diagrams 

A(B°^K-K + ) = P' + T' + ^P^ + P 1 ^ 

A(B Q ^K°7Z Q ) = -±(P'-C'-P^-^ W + P^, (2) 

A(B-^K°JT) = P'- l -P^+A' + P' A , 

A(B-^K-n°) = + T> + C + ^ w + |i^+A' + f A ), 

where T, C, E, A, Pew and Pg W are color-allowed tree, color-suppressed tree, W-exchange, W-annihilation, 
color-allowed and color-suppressed electroweak penguin amplitudes, respectively, and P A is the penguin- 
induced weak annihilation amplitude. We use unprimed and primed symbols to denote A5 = and |A5| = 1 
transitions. We notice that if C' , P EW and A' are negligible compared with T', it is clear from Eq. © that the 



We have included chirally enhanced but power suppressed penguin contributions. Numerically, they are of order 

l/ml 
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decay amplitudes of K 71° and K n + will be the same apart from a trivial factor of 1 / V2. Hence, one will 
expect that A C p(K~K°) « A C p(K~ while they differ by 5.3a experimentally, AA Kjz = Acp{K~% q ) - 
A C p{K-n + ) = 0.148 ± 0.028 We also notice that the decay B -> K r] has a world average -0.37 ± 
0.09 for A C p(K~r\) flflB] different from zero by 4.1 standard deviations. 

Since in the heavy quark limit, CP asymmetries of the K~7l°, K 7], 7T 77 , 7r°7r° modes have the correct 
signs when compared with experiment, the B-CP puzzles mentioned here are relevant to QCDF and may 
not occur in other approaches such as pQCD. In this work, we shall show that soft power corrections to the 
color-suppressed tree amplitude will bring the signs of A C p back to the right track. As a bonus, the rates of 
B° — ► n°n°, p°n° can be accommodated. 

2. The aforementioned direct CP puzzles indicate that it is necessary to consider subleading power cor- 
rections other than penguin annihilation. For example, the large power corrections due to P' A cannot explain 
the AAk k puzzle as they contribute equally to both B — ► K~% and B° — > K~n + . The additional power 
collection should have little impact on the decay rates of penguin-dominated decays but will manifest in the 
measurement of direct CP asymmetries. Note that all the "problematic" modes receive a contribution from 
c O = (f> + P^ff. Since A(B~ -> K~K°) « t' + d + p' and A(B Q -> K~% + ) °c t' +pf with t' = T' + P E C W and 
p' = P' — \P^w + P' A , we can consider this puzzle resolved, provided that d ft' is of order 1.3 ~ 1.4 with a 
large negative phase (naively \d /t'\ ~ 0.9). There are several possibilities for a large d: either a large color 
suppressed C' or a large electroweak penguin or a combination of them. Various scenarios for accom- 
modating large C fl @, 0, [ill, H S Q, 0] or 0,0] have been proposed. To get a large C', one 



can appeal to spectator scattering or final-state rescattering (see discussions below). However, the general 
consensus for a large P^ w is that one needs New Physics beyond the Standard Model. In principle, one can- 
not tell the difference of these two possibilities in penguin-dominated decays as it is always the combination 

c' = C' + 

^ew tnat enters i nto tne decay amplitude except for the decays involving r\ and/or X]' in the final 



state where both d and P^ w present in the amplitudes 11811 . Nevertheless, the two scenarios can lead to very 



distinct predictions for tree-dominated decays where Pew <C C as the electroweak penguin amplitude here 
does not get a CKM enhancement. The decay rates of 5° — > n ^ ',p 7T will be substantially enhanced for 
a large C but remain intact for a large P EW . Since Pew <C C in tree-dominated channels, CP puzzles with 
7t~ri and % 7t cannot be resolved with a large P EW . Therefore, it is most likely that the color-suppressed 
tree amplitude is large and complex. Motivated by the above observation, in this work we shall consider 
the possibility of a large complex c?2, the parameter for describing the color-suppressed tree topology, and 



parameterize power corrections to a.2, as 2 



fla-of^l+pce*), (3) 



with the unknown parameters Or and <pc to be inferred from experiment. 



The reader is referred to j22ll for details. We shall first consider soft corrections to weak annihilation 
dictated by the parameters and <j>A- A fit to the data of two-body hadronic decays of B Q and B mesons 



2 We use NLO results for «2 in Eq. (0 as a benchmark to define power corrections. The NNLO calculations of 
spectator-scattering tree amplitudes and vertex corrections at order a 2 have been carried out in 11911 and 12011 . 



respectively. While NNLO corrections can in principle push the magnitude of 02(^7^) up to the order of 0.50 by 



lowering the value of the Z?-meson parameter Xb, the strong phase of «2 relative to a\ cannot be larger than 15° ||21 ] 
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within QCDF yields the values 



pi rj 1.10, 1.07, 0.87, 



$ » -50°, -70°, -30°, (4) 

for B — > PP, Py, VP respectively, where the superscript "0" of p A and (j) A indicates that they are the default 
values we shall use in this work. Basically, this is very similar to the "scenario S4" presented in [4|]. For 
the annihilation diagram we use the convention that Mi (M2) contains an antiquark (a quark) from the weak 
vertex. Since the penguin annihilation effects are different for M\=P and M\ = V, the parameters p A and 
<f) A are thus different for B — > PV and B — ► VP. 

Branching fractions and direct CP asymmetries for some selective B — ► PP decays are shown in Table U 
The theoretical errors correspond to the uncertainties due to variation of (i) the Gegenbauer moments, the 
decay constants, (ii) the heavy-to-light form factors and the strange quark mass, and (iii) the wave function 
of the B meson characterized by the parameter Ag, the power collections due to weak annihilation and hard 
spectator interactions described by the parameters Pa.h, $a,h, respectively. To obtain the errors shown in 
Table H we first scan randomly the points in the allowed ranges of the above nine parameters (specifically, 
the ranges p° A - 0.1 < p A < p? ±0.1, $ - 20° < (j) A < 0^ + 20°, < p H < 1 and < <fc < 2n are used in 
this work) and then add errors in quadrature. More specifically, the second error in the table is referred to the 
uncertainties caused by the variation of p A n and 0a,//, where all other uncertainties are lumped into the first 
error. Power corrections beyond the heavy quark limit generally give the major theoretical uncertainties. 

For p c ~ 1.3 and C —70°, we find that all the CP puzzles in B — > PP decays are resolved as shown 
in fourth column of Table J] The corresponding a^'s are 

a 2 {nn) « 0.60e~'' 55 °, a 2 {Kn) « 0.51^'' 58 °. (5) 

They are consistent with the phenomenological determination of C^ /T^ ~ az/ai from a global fit to the 
available data [18]. Due to the interference between the penguin and the enhanced color-suppressed am- 
plitudes with a sizable strong phase, it is clear from Table U that theoretical predictions for Aqp now agree 
with experiment in sign even for those modes with the measured Acp less than 3 a in significance. As first 
emphasized by Lunghi and Soni |23|], in the QCDF analysis of the quantity AAkk, although the theoreti- 
cal uncertainties due to power corrections from penguin annihilation are large for individual asymmetries 
Acp(K~n Q ) and Acp(K~n + ), they essentially cancel out in their difference, rendering the theoretical pre- 
diction more reliable. We find AA Kn = (123^1111))%, while it is only (1.9+04+} q)% in the absence of 
power corrections to the topological amplitude "C" or 02- 

For the direct CP asymmetry of B° -> K°n°, we predict A C p(K°n°) = (-10.6+3 7+^)%. Experimen- 
tally, the current world average —0.01 ±0.10 is consistent with no CP violation because the BaBar and Belle 
measurements, — 0.13±0.13±0.03 [24] and 0.14 ±0.13 ±0.06 02511 respectively, are opposite in sign. Nev- 
ertheless, there exist several model-independent determinations of this asymmetry: one is the SU(3) relation 
Ar(7r°7r°) = —Ar(K°7l°) lE^ . and the other is the approximate sum rule for CP rate asymmetries |27] 

Ar(^^ + )±Ar(^ ^") «2[Ar(/r7r )±Ar(£ 7r°)], (6) 

based on isospin symmetry, where Ar(^^) = F(B — > Kn) — Y{B — > Kit). This sum rule allows us to extract 
Acp(K Q n°) in terms of the other three asymmetries in K~7t + ,K~7i Q ,K K~ modes that have been mea- 
sured. From the current data of branching fractions and CP asymmetries, the above SU(3) relation and 
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TABLE I: CP-averaged branching fractions (in units of 10~ 6 ) and direct CP asymmetries (in %) of some 
selective B — > PP decays obtained in QCD factorization for three distinct cases: (i) without any power 
corrections, (ii) with power corrections from penguin annihilation, and (iii) with power corrections to both 
penguin annihilation and color-suppressed tree amplitudes. The parameters Pa and 0a are taken from Eq. 
dU, pc = 1.3 and <p c = —70°. Sources of theoretical uncertainties are discussed in the text. 



Modes 



W/o p A ,c,<l>A,c 



With p A , (j) A 



With p A C , a ,c 



Expt. [6J 







13 1+5-8+0.7 
^ -3.5-0.7 


19 3+ 7 - 9 +F 

^"-'-4.8-6.2 


in 0+7.9+8.2 


19.4 ±0.6 


^fB° — 

Ck> id 


* iY JL J 


<r <r+2.8+0.3 
3 •-'-1.7-0.3 


o 4+3.8+3.8 
°- H -2.3-2.9 


^+3.8+3.8 
°- u -2.2-2.9 


q 4. f. 


J?(£~ _ 


+ K°x~) 


14 q+6.9+0.9 


91 7+9.2+9.0 
Z1 - '-6.0-6.9 


91 7+9.2+9.0 
Z1 -' -6.0-6.9 


23.1 ±1.0 


£%{B~ - 


+ K-K ) 


q 1+3.6+0.5 
-2.3-0.5 


lz -°-3.0-3.7 


1 9 <r+4.7+4.9 
iz - J -3.0-3.8 


12.9±0.6 


S8{B~ - 


* K~n) 


1 -"-0.7-0.4 


2 4+1,+H 

^•^-1. 1-1.0 


-^-l. 1-1.0 


2.3 ±0.3 a 


ggfif — 


► 7T+7T~ ) 


u -0.6-0.4 


7 a+0.4+0.7 
' - u -0.7-0.7 


7 A+0.4+0.7 
' - u -0.7-0.7 


5 16 ±0 22 


m(t - 


► 7r 7T ) 


n 49+0.29+0.18 
U-H -0.1 1-0.08 


a ^9+0.26+0.21 
u - JZ -0. 10-0.10 


, 1+1.0+0.7 
i- 1 -0.4-0.3 


1.55 ±0.19* 


8(Br - 


■> 7l~7l°) 


4 q+0.9+0.6 
^-0.5-0.3 


4 q+0.9+0.6 
^-0.5-0.3 


<r 0+2.2+1.4 
J - y -\. 1-1.1 


5 5Q+0.41 
~'-~' ; '-0.40 


8{Br - 


+ K-ri) 


A 4+0.6+0.4 
H - H -0.3-0.2 


4 r+0.6+0.5 
H - J -0.3-0.3 


<r a+1.2+0.9 
J - u -0.6-0.7 


4.1±0.3 fl 






4 n+0.6+1.1 
H - u -0.7-l.l 


7 /.+1.7+4.3 
'•^-1.5-4.8 


7 4+1.7+4.3 
'•^-1.5-4.8 


-9-8+1;? 


A C p(B° 


-Z°JE°) 


_ 4 a+1.2+3.5 
^• u -l. 8-3.0 


a 7^+1.88+2.56 
u - /J -0.94-3.32 


_i 0fi +2.7+5.6 
1U ' D -3.8-4.3 


— 1 ± 10 


A CP (B- 




a 79+O.O6+O.O5 

u - ' z -o.os-o.05 


a 9O+0.03+0.09 
u - zo -0.03-0.10 


n 90+0.03+0.09 

u - zo -0.03-0.10 


0.9 ±2.5 


A CP (B- 




7 o + t.6+2.3 
'•-'-1.2-2.7 


c+1.3+4.9 

J "'-1.8-4.6 


4 Q+3.9+4.4 
^•^-2.1-5.4 


5.0±2.5 


A C p(B- 


-+K-TI) 


99 ,+7.7+14.0 
zz - i -16.7- 7.3 


19 7+7.7+13.4 
1Z " '-5.0-15.0 


1 1 n + 8-4+14.9 
il - u -21.6-10.1 


-37±9 fl 


A C p{B° 


->• 71+71") 


£ 9+0.4+2.0 
°- -0.5-1.8 


17 a+1.3+4.3 
—1.2—8.7 


17 Q+1.3+4.3 
1 '- u -1.2-8.7 


38±6 


A C p(B° 


«°jf°) 


oo 4+ 6.8+34.8 
10.6-37.7 


Of, n+8.4+48.5 
^ v - y -6.0-37.5 


cn 9 + 14.8+30.3 
D 1 - z -20.8-34.6 


43+25 
«_24 


A C p(B- 


— >■ 7r~7r°) 


u - uo -0.01-0.02 


-0 06+ 01 + 01 
u - uo -0.01-0.02 


n 1 1+O.OI+O.O6 
u - 11 -0.01-0.03 


6±5 


A C p(B- 


->JT*l) 


_i , 4+1.1+2.3 
A1 - H -1.0-2.7 


1 1 4+0.9+4.5 
A1, -0.9-9.1 


A+2.4+ 8.4 
J - u -3.4-10.3 


-13±7 fl 



"We have taken into account the new measurement oLE — ► (K ,n )r\ (7j to update .the average. 
6 This is the average of 1.83 ±0.21 ±0.1 3 by BaBar J3J] and 1.1 ±0.3 ±0.1 by Belle J32J. If an S factor is included, 



the average will become 1 .55 ± 0.35 . 



CP-asymmetry sum rule lead to A CP (^ 7r ) = -0.073 and A C p(K°n ) = -0.15 ±0.04, respectively. 
An analysis based on the topological quark diagrams also yields a similar result —0.08 ~ —0.12 J28[]. All 
these indicate that the direct CP violation Acp(K°n°) should be negative and has a magnitude of order 
0.10. As for the mixing-induced asymmetry S„p K , it is found to be enhanced from 0.76 to 0.791°) 04^0 04 
when pc and (j>c sue, turned on, while experimentally it is 0.57 ±0.17 J6 . The discrepancy between theory 
and experiment for S n o Ks is one of possible hints of New Physics (29J]. Our result for S K o Ks is consis- 
tent with II ll Il2l 1 1 3TI where soft corrections to ci2 were considered, but not with 11411 where S^k - ~ 0.63 
was obtained. A correlation between S^og and Acp{ft Ks) nas been investigated recently in 13011 . For the 
mixing-induced asymmetry in B — > 7T+7T - , we find S n + n - = — 0.69^° io^oo9' i n accordance with the world 
average of -0.65 ±0.07 fl]. 

From TableUwe see that power corrections to the color-suppressed tree amplitude have almost no impact 
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on the decay rates of penguin-dominated decays, but will enhance the color-suppressed tree dominated decay 
B — ► n° 71° substantially owing to the enhancement o f \ao \ ~ ^(0.6). Notice that the central values of the 
branching fractions of this mode measured by BaBar [31] and Belle 13211 are somewhat different as noticed 
in Table HI It is generally believed that direct CP violation of B — > k~iz is very small. This is because the 
isospin of the n~ 71° state is / = 2 and hence it does not receive QCD penguin contributions and receives only 
the loop contributions from electroweak penguins. Since this decay is tree dominated, SM predicts an almost 
null CP asymmetry, of order 10~ 3 ~ 10~ 4 . What will happen if ci2 has a large magnitude and strong phase ? 
We find that soft corrections to the color-suppressed tree amplitude will enhance Acp(k~x q ) substantially to 
the level of 2%. Similar conclusions were also obtained by the analysis based on the diagrammatic approach 
|18|]. However, one must be very cautious about this. The point is that power corrections will affect not only 
ai, but also other parameters a,- with i / 2. Since the isospin of n~7i° is / = 2, soft corrections to aj and 
a,- must be conspired in such a way that 7i~ 71° is always an / = 2 state. As explained below, there are two 
possible sources of power corrections to aj. spectator scattering and final-state interactions. For final-state 
rescattering, it is found in J33I1 that effects of FSIs on Acp(k~k q ) are small, consistent with the requirement 
followed from the CPT theorem. In the specific residual scattering model considered by one of us (CKC) 
lllin . 7l~7l° can only rescatter into itself, and as a consequence, direct CP violation will not receive any 
contribution from final-state interactions. Likewise, if large pn and fa are turned on to mimic Eq. (O, we 
find Acp(7l~7l°) is at most of order 10~ 3 . This is because spectator scattering contributes to not only a 2 but 
also a\ and the electroweak penguin parameters a 7 _ 10 . Therefore, a measurement of direct CP violation in 
B — ► 7r~7T still provides a nice test of the Standard Model and New Physics. 

In order to explain CP violation in the decay B — > K v\ , we shall elaborate it in more detail. Its decay 
amplitude is given by |Q] 



V2A(B~ —>K~rj) = A 



+ V2A 



8 pu OC2 + 20;^ + — GS^EW 



3.EW 



(7) 



+ V2A Kllc [5 pc a 2 + <] +A TJqg \s pu (ai + j3 2 ) + a p A + j8 3 p + /3 3 P EW 



where the flavor states of the 77 meson, 
respectively. The reader is referred to 



q = (uii + dd) I s/2, ss and cc are labelled by the r\ q , r\ s and T]°, 
for other notations. The physical states 7], tj' and r\ c can be 
expressed in terms of flavor states r\ q , t] s and T]°. Since the two penguin processes b — > sss and b — ► sqq 
contribute destructively to B — > Kr\ , the penguin amplitude is comparable in magnitude to the tree amplitude 
induced from b — > usu, contrary to the decay B — ► Kr]' which is dominated by large penguin amplitudes. 
Consequently, a sizable direct CP asymmetry is expected in B — > K r\ but not in K r\' IB411 . 

Quantities relevant to the calculation are the decay constants and denned by (0\qj^Y5q\ri) = 

ifq/v2qn, (0\s YiiY^s \ri) = if^q^ and (0|c7^75c|t]) = if^q^, respectively. A straightforward perturbative 
calculation gives 13511 

f n = -T^ Ij k- (8) 



12m2 y/2 

For the decay constants f% and we shall use the values = 107 MeV and 
B6I1 with the convention of f n = 132 MeV. Although the decay constant : 



-1 12 MeV obtained in 
2 MeV is much smaller than 
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fq ,s , its effect is CKM enhanced by V c bV* s /(V u bV* s ). In the absence of power corrections to a 2 ,A C p(K~r\) is 
found to be 0.127 (see Table|I]). When pc and <pc are turned on, Acp(K~r\) will be reduced to 0.004 if there 
is no intrinsic charm content in the tj. When the effect of /£ is taken into account, Acp(K~r]) finally reaches 
at the level of —11% and has a sign in agreement with experiment. Hence, CP violation in B — > K r\ is 
the place where the charm content of the 7] plays a role. 

We add a remark here that the pQCD prediction for Acp{K~r\) is very sensitive to m qq , the mass of 



the r\ q , which is generally taken to be of order m n . It was found in [37] that for m qq = 0.14, 0.18 and 
0.22 GeV, Acp{K~r]) becomes 0.0562, 0.0588 and —0.3064, respectively. There are two issues here: (i) 
Is it natural to have a large value of m qq ? and (ii) The fact that Acp(K^rj) is so sensitive to m qq implies 



that the pQCD prediction is not stable. Within the framework of pQCD, the authors of 08H rely on the 
NLO corrections to get a negative CP asymmetry and to avoid the aforementioned issues. At the lowest 
order, pQCD predicts Acp(K~i]) « 9.3%. Then NLO corrections will change the sign and give rise to 



A C p{K-i)) = (-n.it?&)%tm 



As for the decay B — > 7177, it is interesting to see that penguin annihilation will flip the sign of 
Acp(7i~ri) into a wrong one without affecting its magnitude (see Table JJ). Again, soft corrections to a 2 
will bring the CP asymmetry back to the right track. Contrary to the previous case, the charm content of 
the r\ here does not play a role as it does not get a CKM enhancement relative to the non-charm content 
of the r\. Our result of Acp(n~r]) = — 0.05 jig 1? is consistent with the measurement of —0.13 ±0.07. For 



comparison, the pQCD approach predicts -0.37^ ^ H39Q and SCET gives two solutions 140], 0.05 ±0.29 
and 0.37 ± 0.29 with signs opposite to the data. 

3. What is the origin of power corrections to a 2 ? There are two possible sources: spectator scattering 
and final-state interactions. The flavor operators af are basically the Wilson coefficients in conjunction 
with short-distance nonfactorizable corrections such as vertex corrections, penguin contractions and hard 
spectator interactions. In general, they have the expression 0, 01 

afiMM) = ^ i . + ^^ i .(M 2 )±^^[yKM 2 )±^// ; (M 1 M 2 )]±^(M 2 ), (9) 

where i = 1 , - • ,10, the upper (lower) signs apply when i is odd (even), a are the Wilson coefficients, 
C F = (N% - l)/(2N c ) with N c = 3, Ni(M 2 ) = for i = 6,8 and equals to 1 otherwise, M 2 is the emitted 
meson and Mi shares the same spectator quark with the B meson. The quantities Vj(M 2 ) account for vertex 
corrections, Hi(M\M%) for hard spectator interactions with a hard gluon exchange between the emitted 
meson and the spectator quark of the B meson and P, (M 2 ) for penguin contractions. A typical hard spectator 
term Hi(M\M 2 ) has the expressions |Q,0]: 

Hi ( Ml M 2 ) = IMmJm ; ? f'dxdy ( W + rf , (10 ) 

for i = 1-4,9,10, where X {m ^ M ^ is the factorizable amplitude for B — > M\M 2 , x = 1 -x, X B is the fist 
inverse moment of the B meson light-cone wave function and 

= . v 2m i v - v r» = 2n ;\ f *y . (id 

Power coiTections from the twist-3 amplitude are divergent and can be parameterized as 

-1 dy mp 



X H = I Z^=ln^i(l+p^«), (12) 
y A ft 
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II 




FIG. 1: Contribution to the color-suppressed tree amplitude of B — > K~n from the weak decay B — > 
^"17' followed by the final-state rescattering of if" 77' into K°n Q . This has the same topology as the color- 
suppressed tree diagram. 



Since c\ ~ i^(l) and C9 ~ ^(—1.3) in units of a em , it turns out that spectator scattering contributions to a,- 
are usually small except for a 2 and a\$ which are essentially governed by hard spectator interactions 14111 . 
The value a 2 (Kn) m 0.51 e~' 5& ° corresponds to p# « 4.9 and <pH ~ —77°. 3 Therefore, there is no reason to 
restrict p# to the range < p# < 1. 

A sizable color-suppressed tree amplitude also can be induced via color-allowed decay B — > K~r\' 
followed by the rescattering of K r\' into K~n° as depicted in Fig. 1. Recall that among the 2-body B 
decays, B — ► Kr\' has the largest branching fraction, of order 70 x 10~ 6 . This final-state rescattering has 
the same topology as the color-suppressed tree diagram 13311 . One of us (CKC) has studied the FSI effects 
through residual rescattering among PP states and resolved the B-CP puzzles jllH . 



4. Power corrections to a 2 for B — ► VP and B — > VV are not the same as that for B — ► PP as described 
by Eq. ([5]). From Table [II] we see that an enhancement of a 2 is needed to improve the rates of B — > p°7r° 
and the direct CP asymmetry of B° — > K* 0r q. However, it is constrained by the measured rates of p°7T~ and 
p~7T° modes. This means that pc(V-P) is preferred to be smaller than pc(PP) = 1.3 . In Table HI] we show 
the branching fractions and CP asymmetries in B — ► VP decays for pc(VP) = 0.8 and <j> c (VA) = —80°. The 
corresponding values of a 2 {VP) are 



a 2 (np) «0.40e 
a 2 {pK) «0.36e 



-(52° 



a 2 {pn) 
a 2 (%K* 



0.38e 
%0.39e 



-i52° 



-i'51° 



(13) 



It is clear from Table HT1 that in the heavy quark limit, the predicted rates for B —> K*n are too small by 
a factor of 2 ~ 3, while SS(B — > ^p) are too small by (15 ~ 50)% compared with experiment. The rate 
deficit for penguin-dominated decays can be accounted by the subleading power corrections from penguin 
annihilation. Soft corrections to a 2 will enhance 3§(B — ► p°7i°) to the order of 1.3 x 10 -6 , while the BaBar 
and Belle results, (1.4 ±0.6 ±0.3) x 10~ 6 H and (3.0 ± 0.5 ± 0.7) x 10~ 6 H respectively, differ in 
their central values by a factor of 2. Improved measurements are certainly needed for this decay mode. As 



3 As pointed out in I2lll42ll . a smaller value of Xb of order 200 MeV can enhance the hard spectator interaction [see 
Eq. ( fTOl il and hence a 2 substantially. However, the recent BaBar data onB^ yiv M43I1 seems to imply a larger Xb 
(> 300 MeV at the 90% CL). In this work we reply on pc and <j) c to get a large complex a 2 . 
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TABLE II: Same as Table U except for some selective B — > VP decays with p c = 0.8 and <pc = —80°. 



iviooes 




w/o Pa,Ci vm,c 


witn Pa, 


wiua Pa,c ; 9a ,c 


u vnr r£n 


ix> Id — 


►A p ) 


£ ir+5.4+0.4 
°- D -2.6-0A 


£+5.7+7.4 

°- D -2.8-4.5 


^+5.7+7.4 

°- D -2.8-4.5 


— 1.1 


tx>\p — 




4 7+3.3+0.3 
^-'-1.7-0.3 


5 r+3.5+4.3 
J - J -l. 8-2.8 


5 4+3.3+4.3 
J -^-1.7-2.8 


5 4+0.9 
5 - 4 -1.0 


(D — 


+ A p ) 


^ 5+6.1+0.7 
J - J -2.8-0.5 


7 0+6.3+7.3 
'•°-2.9-4.4 


7 0+6.3+7.3 
'•°-2.9-4.4 


Q A+1.5 


(31(R— 


+ A p ) 


, q+2.5+0.3 
1 -^-1.0-0.2 


o o+2.6+2.9 
-'•-'-1.1-1.7 


o 5+2.9+2.9 
^ "'-1.2-1.8 


o ,+0.48 
- 5 - 8i -0.46 


tx>\p — 


►A K ) 


o 7 +0.5+0.4 
-'■ ' -0.5-0.4 


9 + 1.0+3.7 
^• z -1.0-3.3 


a 9+1.0+3.7 

^• z -1.0-3.3 


IU.j ±1.1 


ix>[B — 


►A 71 ) 


i ,+0.2+0.2 
1 - 1 -0.2-0.2 


o 5+0.4+1.7 
-'• J -0.5-1.5 


o 5+O.4+I.6 
—0.4— 1 .4 


a 1 n n 
LA ±0. / 


tx>\p — 


* A K ) 


A n+0.7+0.6 
^■ u -0.9-0.6 


m 4+1.3+4.3 

1U - 4 -1.5-3.9 


m 4+1.3+4.3 
1U - 4 1.5 3.9 


Q Q+0.8 


(3>(Tt— 


. TS* -77-0 \ 

■> A % ) 


t 9+0.4+0.3 
-••^-0.4-0.3 


£ O+0.7+2.3 
D -°-0.7-2.2 


^ 7+0.7+2.4 
D - '-0.7-2.2 


u.y ± z.j 


OS Yr° 


Y*®m\ 

►A T]j 


i i ft +6.9+1.7 
1LU -3.5-1.0 


1 5 4+7.7+9.4 
15.4_ 40 _ 71 


1 5 £+7.9+9.4 

^•— 4.1— 7.1 


ij.y ± i.u 


c£> Id — 


► n°ir 0N l 

► P K ) 


7fi+0.96+0.66 
u - /D -0.37-0.31 


n 50+0.88+0.6O 
u -' o -0.32-0.22 


, 0+I.7+I.2 
1 "'-0.6-0.6 


Z.U It U.J 


OfilR- 


* p 7Z ) 


11 fi +1.2+0.9 
ii - u -0.9-0.5 


11 O + 1.3+1.0 
11 -°-0.9-0.6 


1 . + 1.8+1.4 
1 18 1.1 1.4 


iUy -1.5 


OfllJ)- 

C&\D — 


*P % ) 


o 9 + 1.8+1.2 
°- z -0.9-0.6 


5+I.8+I.2 
°- J -0.9-0.6 


7+2.7+1.7 
8 - '-1.3-1.4 


o+1.2 
3 


^> \D — 


' P 71 ) 


ic o + LO+0.5 
1J --'-1.5-0.9 


15 q+1. 1+0.9 
13y 1.5-1.1 


15 q+1. 1+0.9 


1^7 + 1 S 
1 J . / 31 l.O 


O#(R — 

Cxy ID 


P 71 1 


o 4 +0.4+0.3 
"•^-0.7-0.5 


a 9+0.4+0.5 
*- z -0.7-0.7 


n 9+0.4+0.5 

^• z -0.7-0.7 


7^+19 


A fR° 


_v Y n + \ 
' K P ) 


i O+0.7+3.8 
i --'-0.3-3.8 


11 Q+ll.5+19.6 
■ 51 - y -11.0-12.7 


il Q+ll.5+19.6 
- ,i -^-11.0-12.7 


1^4-6 
1 J ± U 


A (R° 


-^A p j 


^ o + l. 1+4.9 
b - 8 -1.2-4.9 


5 n+3.2+6.0 
J - u -6.4-4.5 


7+1.2+8.7 
°- -1.2-6.8 


6 4-10 
± iz 


A ( R 


—►A p J 


n 9 4 +0. 12+0.08 
"■^-0. 15-0.07 


n 97+0.19+0.46 
-0.27-0.17 


n 97+0.19+0.46 
u - z/ -0.27-0.17 


19 4-17 
— IZ ± 1 / 


a „„(r— 


i Y-n®\ 

— >A p J 


o O+3.5+7.0 
°- J -0.9-7.0 


5^ 5+16.1+30.0 
JD - J -18.2-22.8 


45 4+17.8+31.4 
19.4-23.2 


44 -17 


A (R° 


— > A 71 J 


1 ^ A+0.9+4.5 
l:) - D -0.7-4.7 


19 1+0.5+12.6 
-0.5-16.0 


19 1+0.5+12.6 
lz - 1 -0.5-16.0 


—23 ± 


a m° 

A C p(B 


v r>*0_0\ 
—>■ A 7T J 


19 n+2-4+11.3 

-12.0_ 4 6 _ 7 6 


n O7+1.71+6.04 
u -°' -0.89-6.79 


m 7+1-8+9.1 
1U - '-2.8-6.3 


— 1 J ± 11 


A ( D 

Acp{B 


— >• A 7T j 


n Q7+0.11+0.12 
v y '-0.07-0.il 


n OQ+0.04+0.10 
u '-'^-0.03-0.12 


O oq+0.04+0.10 
u -' y -0.03-0.12 


3.2 ±0.4 


A CP [B 


— ► A K ) 


17 r+2.0+6.3 
1 ' - J 1.3-8.0 


^ 7+O.7+II.8 
D - '-1.1-14.0 


1 £+3.1 + 11.1 

1 -°-1.7-14.3 


4 ± zy 


A C p(B° 




t A+0.4+1.9 
-'• u -0.4-1.8 


n 90+0.51+2.00 

u - z,u - 1.00- 1.21 


o 5+0.4+2.7 
-'• J -0.5-2.4 


19±5 


A C p(B° 


-p°7r°) 


9 ^,+2.4+9.9 
' 6 -' 3 -3.7-9.2 


ti 5+13.3+21.5 
- ,1 - J -12.5-30.9 


1 1 n+5.0+23.5 
li - u -5.7-28.8 


-30 ±38* 


A C p(B- 


^p-7T°) 


5 A +0.4+2.0 
J - H -0.3-2.1 


i 6 3+I.I+ 7.1 
iD-J_1.2-l0.5 


Q 7+2.I+ 8.0 
^- ' 3.1-10.3 


2± 11 


A C p(B- 


^p°7T-) 


^ 7+0.5+3.5 
u - -0.8-3.1 


1Q s+1.7+12.6 
i ^-°-1.2-8.8 


Q n+3.4+11.4 

y -°-2.6-10.2 


18+ 7 


A C p(B° 


-^p-% + ) 


o ^+0.2+1.0 
-'• J -0.2-0.9 


4 4+0.3+5.8 
^-0.3-6.8 


4 4+0.3+5.8 
^•^-0.3-6.8 


11 ±6 


A C p{B° 


->p + 7T) 


6 +ai+2 - 2 

u ' u -0.1-2.2 


99 7+0.9+8.2 
zz -'-l. 1-4.4 


99 7+0.9+8.2 
zz -'-l. 1-4.4 


— 18± 12 



! 'If an S factor is included, the average will become 2.0 ± 0.8. 
b This is the average of 10 ± 40 ± 53 by BaBar 144D and -49 ± 36 ± 28 by Belle 145J 



for direct CP asymmetries, we see that penguin annihilation will flip the sign of Acp(K~ p°) into the right 
direction. Power corrections to the color-suppressed tree amplitude are needed to improve the prediction for 
Acp(K*°T]). Our prediction is of order 0.035 to be compared with the experimental value of 0.19 ±0.05. 

The pQCD prediction of A C p(K*°T]) ~ 0.0057 H is too small, while the SECT result of 0.01 H 

has a wrong sign. For Acp(K°p°), it gets a sign flip after including soft effects on a%. Our prediction is 
(8.7+jJ)%, while it is 0.06 ±0.20 experimentally. Defining AA K * n =A C p(K*-71°) -A C p{K*-% + ) in analog 
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to AA Kn , we predict that AA K * n = (13.7+j 4+g|)%, while it is naively expected that K*~K° and K*~K + 
have similar CP-violating effects. It is of importance to measure CP asymmetries of these two modes to 
test our prediction. For mixing-induced CP violation, we obtain AS$k s = , AS m K s = 0- 08 



and AS p o Ks = — 0.17+Jyjg j22ll . where ASf = —i\fSf — sin2/3. It turns out that soft corrections to c?2 have 
significant effects on the last two quantities. 

As for B — > VV decays, we notice that the calculated fi° -> p°p° rate in QCDF is ^(B° ->• p°p°) 



x 10" 6 for p c = 150|1, while the world average is (0.73t££g) x 1CT 6 fej. Therefore, soft 



+0.27 \ v 1n -6 



power correction to ci2 or PciVV) should be small for B° — > p°p°. Consequently, a pattern follows: Effects 
of power coiTections on «2 are large for PP modes, moderate for VP ones and very small for VV cases. 4 



This is consistent with the observation made in H 1 3TI that soft power correction dominance is much larger for 
PP than VP and VV final states. It has been argued that this has to do with the special nature of the pion 
which is a qq bound state on the one hand and a nearly massless Nambu-Goldstone boson on the other hand 



1 13J]. The two seemingly distinct pictures of the pion can be reconciled by considering a soft cloud of higher 
Fock states surrounding the bound valence quarks. From the FSI point of view, since B — ► p + p~ has a rate 
much larger than B — > 7l + 7i~, it is natural to expect that B — ► 71° 71° receives a large enhancement from the 
weak decay B — > p + p~ followed by the rescattering of p + p~ to k°k° through the exchange of the p particle. 
Likewise, it is anticipated that B — > p°p° will receive a large enhancement via isospin final-state interactions 
from B — > p + p . The fact that the branching fraction of this mode is rather small and is consistent with the 
theory prediction implies that the isospin phase difference of 8q and 8% and the final-state interaction must 



be negligible 115 lh - 



5. B-CP puzzles arise in the framework of QCD factorization because power corrections due to pen- 
guin annihilation, that account for the observed rates of penguin-dominated two-body decays of B mesons 
and direct CP asymmetries A C p(K~n + ), A C p(K*~ % + ), A C p(K~p°) and A C p(n + K~), will flip the signs 
of direct CP-violating effects in B — > K~n°,B~ — > K 7],B — > K~r\ and B° — > n°n Q to wrong ones 
when confronted with experiment. We have shown that power corrections to the color-suppressed tree 
amplitude due to hard spectator interactions and/or final-state interactions will yield correct signs again 
for aforementioned CP asymmetries and accommodate the observed n n () and p°7r° rates simultaneously. 
CP- violating asymmetries of B — > K r\ can be understood as a consequence of soft collections to 02- 
Acp(K Q 7l°) is predicted to be of order —0.10, in agreement with that inferred from the CP-asymmetry 
sum rule, or SU(3) relation or the diagrammatical approach. For direct CP violation in B — > K*~r],K~r\, 
our predictions are in better agreement with experiment than pQCD and SCET. For B Q — ► K°p°, we ob- 
tained A CP (K°p ) = 0.087+0069- We argued that the smallness of CP asymmetry of B -> n^n Q is not 
affected by the soft coiTections under consideration. For the CP asymmetry difference in K*n modes de- 



4 Since the chiral factor r v x for the vector meson is substantially smaller than r£ for the pseudoscalar meson (typically, 
r X = ^(0-8) and = 0(0.2) at the hard collinear scale jj. = \/Am&), one may argue that Eq. ( fTOb naturally explains 
why the power corrections to 02 is smaller when M\ is a vector meson, provided that soft corrections arise from 
spectator rescattering. Unfortunately, this is not the case. Numerically, we found that, for example, H(K*n) is 
comparable to H(Kn). This is due to the fact that / dxr^<i> m (x) /x is equal to X#r£ for M = P and approximated 
to 3{X H -2)r v x for M = V. 



10 



fined by AAg*„ = A CP (K* 7C°) —A C p(K* n + ), we predict that AA K * n ~ 14%, while these two modes are 
naively expected to have similar direct CP-violating effects. For mixing-induced CP violation, we found 
AS K o Ks = 0A2t° Z , ASt Ks = 0.022^ , AS aKs = 0.17+°$ and AS p o Ks = -0.17+°;°*. 
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